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We have investigated the role of temperature and magnetic effects on the stacking-fault energy 
(SFE) in pure austenitic iron based on density functional theory (DFT) calculations. Using the 
axial next-nearest-neighbor Ising (ANNNI) model, the SFE is expanded in terms of free energies of 
bulk with face-centered cubic (fee), hexagonal close-packed (hep), and double-hep (dhep) structures. 
The free-energy calculations require the lattice constant and the local magnetic moments at various 
temperatures. The earlier is obtained from the available experimental data, while the later is 
calculated by accounting for the thermal magnetic excitations using Monte-Carlo techniques. Our 
results demonstrate a strong dependence of the SFE on the magnetic effects in pure iron. Moreover, 
we found that the SFE increases with temperature. 


I. INTRODUCTION 

Iron and its alloys have played a significant role in the 
development of human civilization. Beside its industrial 
importance, some unique features of iron like its phase 
transitions and magnetic properties have opened inter¬ 
esting fields of research for materials scientist. 

The well-known phase diagram of iron shows that, 
at atmospheric pressure, and at low temperatures, pure 
iron is found in ferromagnetic body-centered cubic (bcc) 
structure^. As temperature rises, at 1043K (770°C), i.e., 
at the Curie temperature of iron, it demonstrates a mag¬ 
netic phase transition from the ferromagnetic to param¬ 
agnetic, while preserving the bcc structure. At 1185 K 
(912°C), iron faces a structural phase transition from the 
bcc to the face-centered cubic (fee) structure. Further 
heating of iron reveals a second structural phase transi¬ 
tion at 1667K (1394°C), through which the fee structure 
changes back to the bcc. Finally, at 1811K (1538°C) 
iron melts^. Although the pure iron is found in the bcc 
structure at room temperatures, the addition of alloying 
elements like manganese and nickel can stabilize its fee 
phase at room temperature'^’^. 

The mechanical properties of steels is influenced by 
their plastic deformations. In fee metals, plastic deforma¬ 
tions may occur through different mechanisms including 
dislocation gliding, twinning (twinning-induced plastic¬ 
ity, TWIP), and phase transformation (transformation- 
induced plasticity, TRIP). The activation of these mech¬ 
anisms has been proven to be governed by the stacking- 
fault energy (SFE)^“®. 

The intrinsic stacking fault (SF) is one of the simplest 
planar defects in the fee crystal lattice'. The fee struc¬ 
ture is formed by the stacking of the close-packed layers in 
the ... ABC ABC ABC... sequence. In this structure, an 
ISF can be considered as the elimination of a close-packed 
layer in the bulk, resulting in the ... AiJCABj ABC... 
sequence, where the removed layer has been denoted by 
a vertical line. This defect is demonstrated in Fig. 1. 
In the neighborhood of the fault, the structure resem¬ 
bles locally the stacking sequence of the hexagonal close- 


packed (hep) structure, as highlighted by bold letters in 
the notation^’® (see Fig. 1). The energy associated with 
a SF, the SFE, is defined as the difference between the 
free energy of a structure with a fault and that of the 
perfect fee structure: SFE = Fsf — Ffec- 



1211J 


FIG. 1. The stacking sequence of an fee structure along its 
[111] direction, with an intrinsic stacking fault at the position 
denoted by the horizontal orange line. The perfect fee stack¬ 
ing order is highlighted in green, while the interruption due 
to the fault is emphasized by the turquoise color. 

There has been numerous experimental works on the 
measuring of the SFE in different austenitic steels®^^^. As 
already discussed by Abbasi^^, Gholizadeh^®, and Reyes- 
Huamantinco^^, the available experimental data for the 
SFE are highly questionable since the reported ranges are 
too broad and have shown a high deviation from average 
amount. Therefore, the development of a theoretical ap¬ 
proach based on quantum mechanics is highly desired. 

In recent years, the SFE and its dependence on differ¬ 
ent parameters have been the subject of many theoretical 
investigations within the framework of the density func¬ 
tional theory (DFT)^®d6^ Some of them have applied 
supercell approaches, which allow for the explicit simu¬ 
lation of the fault as well as for the relaxation of local 
forces, but critically restricts the calculation of magnetic 
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or chemical disorders. For instance, Kibey et have 
studied the dependence of the SFE on the concentration 
of interstitial nitrogen in fee iron. Similarly, Abbasi et 
and Gholizadeh et have studied the influence 
of interstitial carbon of the SFE in fee iron. In these three 
works, the lattice local displacements introduced by the 
interstitial atom and also by the fault are accurately ac¬ 
counted for. On the other hand, using spin-unpolarized 
simulations according to 0 K equilibrium, all magnetic 
interactions in the paramagnetic medium are neglected. 
Abbasi^^ reports only small differences between results 
obtained from tests with nonmagnetic (spin-unpolarized) 
and ferromagnetic calculations, supporting their simpli¬ 
fied nonmagnetic calculations where the influence of the 
magnetic interactions on the qualitative behavior of the 
SFE is assumed negligible. Referring to Abbasi’s tests, 
Gholizadeh^^ states that although a nonmagnetic calcu¬ 
lation may be reliable enough to study the dependence of 
the SFE against the interstitial concentration, calculat¬ 
ing the absolute value of the SFE and also developing a 
complete understanding of the atomic interactions in the 
paramagnetic medium requires accurately accounting for 
the magnetic interactions. 

Other investigations have applied Green’s function 
formalism^®^^^, where the utilization of the disordered 
local moments (DLM)^^ approach and the coherent po¬ 
tential approximation (GPA)^^“^^ allows for the simu¬ 
lation of the magnetic and chemical disorders, respec¬ 
tively. For instance, Vitos et alP'^ studied the tempera¬ 
ture dependence of the SFE in iron-chromium-nickel al¬ 
loys. Vitos concludes that the temperature dependence 
of the SFE is almost totally dictated by the contribu¬ 
tion of the magnetic fluctuations into the free energy. 
Later, Reyes-Huamantinco et al}"^ and Gholizadeh et 
improved Vitos’s approach, particularly by includ¬ 
ing the experimental data for the thermal lattice expan¬ 
sion, and calculated the temperature dependence of the 
SFE in iron-manganese and iron-chromium-nickel alloys, 
respectively. The two works reveal that the temperature 
dependence of the SFE is mainly obtained from the to¬ 
tal energy of the alloy, which is in turn a function of its 
lattice parameter at different temperatures. Therefore, 
in contrast with Vitos’s results, the two works conclude 
that the temperature dependence of the SFE is mainly 
dictated by the lattice thermal expansion. 

Although the published works emphasize that account¬ 
ing for magnetic effects is crucial for understanding the 
phase stability and hence the behavior of the SFE in Fe- 
based alloys, an explicit comparison between quantitative 
results obtained from the paramagnetic calculations and 
those obtained from nonmagnetic (spin-unpolarized) cal¬ 
culations can reveal the magnitude of the magnetic inter¬ 
actions in the SFE. Iron is the dominant element in many 
industrially interested alloys, including those mentioned 
above. Further investigations show that the Fe atoms 
are the main responsible for the magnetic interactions in 
these alloys^^’^^. Moreover, performing calculations for 
pure iron avoids all complexities which are related to an 


alloy compared to an element, like atomic size mismatch, 
local lattice distortion, short-range and long-range or¬ 
ders, Suzuki effect, chemical phase transitions etc. There¬ 
fore, in the current work we compare two sets of the 
SFEs calculated for pure iron in fee phase, one set using 
the methodology introduced by Reyes-Huamantinco et 
and the other set using nonmagnetic calculations. 
Such comparison will answer two main questions: (i) How 
does the SFE change with temperature in the paramag¬ 
netic austenitic iron? and, (ii) How big is the influence 
of the magnetic effects on the SFE in the paramagnetic 
austenitic iron? 


II. METHODOLOGY 

A. The ANNNI Model 




FIG. 2. Primitive cells of three crystal structures used in the 
ANNNI model. Atoms are colored according to their stacking 
position along the [111] direction, (a) depicts the primitive 
cell of the fee structure with only one atomic site. The cu¬ 
bic cell is shown for a better imagination of the lattice, (b) 
shows the primitive cell of the hep structure containing two 
non-equivalent atoms (two atomic sites), (c) represents the 
primitive cell of the dhep strncture with four non-equivalent 
atoms. 


The axial next-nearest-neighbor Ising (ANNNI) model, 
as explained by Cheng et and Denteneer et 
expands the SFE in terms of the free energies of bulk 
unit-cells with different structures. In its second order, 
the ANNNI model results in 

SPE(y) ^ Ehcp(r) + 2j^dhcp(r)_3j,fcc(r) ^ 

.A 

( 1 ) 

where denotes the Helmholtz free energy of a single 
atom in phase (j), A is the area in a close-packed layer 
occupied by a single atomic site, i.e., A = and 

0(3) is the error introduced by neglecting the higher or¬ 
der interactions. 
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B. Temperature Dependence of the Free Energy 

The Helmholtz free energy is defined as 

F(T) = E{T) - TS{T), (2) 

where E{T) and S{T) are the total (internal) energy and 
entropy, respectively. 

The total energy E{T) is calculated using the 
DFT^®d6 j^g temperature dependence originates from 
three sources: (i) the electron distribution over states 
defined by the set of occupation numbers {a^ \ e < ep}, 
(ii) the lattice parameter a(T), and (hi) the average lo¬ 
cal magnetic moment m{T): 

E{T)=E{{a,},a{T),miT)). (3) 

The Mermin functional^^, is applied in finite-temperature 
DFT calculations to account for the temperature depen¬ 
dence of the electron distribution of over states^^. The 
temperature dependence of the lattice parameter is taken 
into account by selecting the lattice parameters according 
to the experimental data for thermal lattice expansions 
(see Fig. 3). Finally, the temperature dependent local 
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FIG. 3. The experimental lattice parameter of pure iron as 
a function of temperature, obtained from high-temperature 
XRD measurements^^. 

magnetic moments are evaluated using a statistical ap¬ 
proach, which is explained in the next subsection. 

C. Longitudinal Spin Fluctuations 

The thermal excitation of the local magnetic moments 
in the paramagnetic DLM state is described using a 
simple model based on the theory of unified itinerant 
magnetism^^. Originally, the theory accounts for both 
transverse and longitudinal spin fluctuations on equal 
footing. However, here it is further simplified so that 
the transverse spin fluctuations, i.e., the fluctuations in 
the orientations of the magnetic moments, always follow 
the completely disordered configuration described in the 
DLM state. The longitudinal spin fluctuations (LSF), 


i.e., the fluctuations in the size of the magnetic moments, 
are obtained by performing a classical Monte-Carlo sim¬ 
ulation over a mapping of the system energetics. These 
system energetics are calculated using the EMTO code. 

For representing the energy of the classical magnetic 
state, Ruban et al.^'^ introduced a magnetic Hamiltonian 
which was later applied to an Fe-Mn alloy by Reyes- 
Huamantinco et In the case of pure iron, the Hamil¬ 
tonian is simplified as 

Hynag.im) = J{m), (4) 

where m is the spatially-averaged local magnetic mo¬ 
ments of iron atoms, and J(rn) is the energy required to 
excite this averaged moment from 0 to the value m in the 
DLM paramagnetic state. After calculating the Hamilto¬ 
nian parameters, J{m), a Metropolis Monte-Carlo tech¬ 
nique is applied to find the temperature-dependence of 
the local magnetic moments^^. 


D. Entropy Contributions 

In a solid, the entropy consists of configurational, vi¬ 
brational, magnetic, and electronic contributions^^ 

S = 5^vib. “t” Sconf. T Smag. T Se\., (5) 

where all terms are associated to a single site in the lat¬ 
tice. 

Currently, there are no available theoretical tools to 
determine the vibrational entropy in paramagnetic ran¬ 
dom alloys. 

The configurational entropy is connected to the disor¬ 
der in the arrangement of atoms of different species in the 
material. Therefore, in a pure element where all atoms 
are equal, it simply vanishes: S'conf. = 0. 

For the ideal paramagnetic state, the magnetic entropy 
is evaluated using the mean-field expression^® 

'S'mag. (T) = kB In {m{T) -b l). (6) 

Here, m(T) denotes the average (over time and space) 
local magnetic moment obtained using the mentioned 
Monte Carlo calculations. 

The electronic entropy can be calculated as®^ 

5ei,(T) = -2A:B !{/(£) In (/(e)) (7) 

+ (1 - /(e)) In (1 - /(e)) }D(e) de, 

where, D{e) and /(e) denote the density of states and the 
finite-temperature Fermi function®®, respectively. The 
Fermi function is a consequence of the Fermi-Dirac 
statistics®®’®® and gives the probability of occupation of 
a state with energy e at temperature T. The Fermi func- 
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tion is defined as 

^ e(e-M)7feBT 

where /i indicates the chemical potential. 


E. The SEE Calculations 

In order to find the total energy as a function of tem¬ 
perature, the temperature-dependent lattice parameter 
and local magnetic moments are used in a set of con¬ 
strained DFT calculations, where the magnetic moment 
is fixed to its corresponding value evaluated using the 
Monte-Carlo method. The free energy is simply found 
by subtracting the entropy contributions from the total 
energy (see Eqs. 2 and 5). The SFE is evaluated using 
these total energies, according to the Eq. 1. 


III. CONCLUSION 

Our results evaluated using both spin-unpolarized and 
spin-polarized calculations show that the SFE increases 
with temperature. Although, spin-polarized calculations 
suggest that the temperature dependence of the SFE van¬ 
ishes at high temperatures. The similar trend of the SFE 
in both calculations shows that the thermal dependence 
of the SFE is mainly originated from the thermal lat¬ 
tice expansion. This is in contrast with the published 
work by Vitos et , but agrees with results found by 
Reyes-Huamantinco^^ and Gholizadeh^ ' . We also find 
that, in order to get the correct phase stability of iron 
and also the correct value of its SFE, one must account 
for magnetic interactions in the system, including the 
spin fluctuations. 
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